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SECTION  I 


INTRODUCTION 

At  the  present  time,  a great  deal  of  attention  is  being  directed  to 
turbine  engine  structural  durability.  Emphasis  is  now  placed  on  either 
extending  or  extracting  the  maximum  useful  life  from  engine  hardware. 
With  major  advances  being  made  in  fracture  mechanics,  the  day  may  soon 
arrive  when  an  engine's  structural  components  will  be  designed  to 
satisfy  a damage  tolerance  criteria  similar  to  that  employed  on  air- 
frames today.  Therefore,  total  understanding  of  the  effect  of  dis- 
crepancies, i.e.,  flaws  and  cracks  on  engine  hardware,  is  necessary 
before  such  an  approach  to  design  can  be  feasible. 

One  area  which  may  be  especially  sensitive  to  these  discrepancies 
is  turbine  engine  blades  (i.e.,  fan,  compressor,  and  turbine).  Engines 
are  designed  to  meet  certain  requirements  which  are  determined  by  blade 
response  (i.e.,  flutter  boundaries,  engine  order  blade  excitations, 
etc.).  Both  blade  flutter  and  frequency  response  are  dependent  upon 
material  and  mechanical  damping  characteristics,  blade  mode  shapes, 
and  natural  frequencies  of  excitation.  This  paper  will  explore  the 
effects  of  cracks  on  the  mode  shapes  and  frequencies  of  blades.  For 
the  purpose  of  portraying  the  relative  effects  of  cracks,  a flat  plate 
with  an  aspect  ratio  of  two  will  serve  as  the  blade  model  for  this 
investigation. 

Much  work  has  been  directed  to  the  solution  of  the  flat  plate 
vibration  problem.  Extensive  documentation  (Reference  l)  of  upper  and 
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lower  bounds  of  solutions  for  various  plate  geometries  has  meticulously- 


defined  what  the  exact  value  of  solution  would  be,  provided  it  could 
be  obtained.  Various  techniques  to  solution  have  been  employed  yield- 
ing reasonable  accuracy,  however,  no  exact  solution  is  available  for  the 
nonuniform  (cracked)  boundary  value  problem.  Probably  one  of  the  most 
powerful  tools  for  the  solution  of  a problem  of  this  nature  is  that  of 
the  finite  element  theory.  The  utility  of  such  an  approach  to  solution 
is  questionable  when  a more  direct,  simplified,  and  economical  means 
of  solution  is  available;  namely,  the  energy  solution.  It  is  unfortu- 
nate that  some  energy  approaches  to  solutions  (i.e.,  Rayleigh,  Galerkin) 
have  masked  the  power  of  this  type  of  approach.  The  major  drawback  of 
these  approaches  has  been  the  need  of  selecting  "appropriate"  shape 
functions  which  satisfy  the  geometrical  constraints  and  approximate 
the  modes  of  vibration.  The  accuracy  of  the  solutions  are  primarily 
dependent  on  these  assumed  functions.  Ritz  (Reference  2)  demonstrated 
that  through  the  use  of  a truncated  series,  properly  employed,  an  upper 
bound  on  the  true  solution  to  deformable  body  problems  could  be  obtained. 

The  benefits  of  using  an  approach  to  solution  as  eluded  to  by  Ritz  can 
best  be  summarized  as  follows  (Reference  3).  "Advantages  of  the  Ritz 
method  lie  in  the  relative  ease  with  which  complex  boundary  conditions 
can  be  handled.  It  is  a powerful  tool  yielding  high  accuracy  in  the 
deflection  analysis  . . . The  Ritz  method  can  be  considered  as  one 
of  the  most  usable  methods  of  higher  analysis  for  solving  complex 
boundary  value  problems  in  the  mathematical  physics."  The  relative 
ease  and  accuracy  of  this  approach  combined  with  the  concepts 

2 

— IU 


associated  with  Hamilton's  Law  of  Varying  Action  provides  an  analytical 
means  of  solution  for  a wide  range  of  applications  (i.e.,  conservative, 
nonconservative,  stationary,  and  nonstationary  motion  of  particles, 
beams,  plates  and  shells  (References  4 through  8). 

This  paper  will  concentrate  on  the  stationary  motion  of  the  simple 
harmonic  vibration  of  thin  flat  plates.  Since  no  in-plane  vibrations 
(membrane)  will  be  considered,  the  out-of-plane  deflection  w can  be 
written  as  a function  of  the  in-plane  coordinates  x,  y,  and  time,  t. 

As  will  be  shown,  the  assumptions  associated  with  linear,  plane  stress 
solutions  may  no  longer  be  valid  to  solve  a highly  nonuniform  vibration 
problem  such  as  in  the  crack  plate  response. 


. 
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SECTION  II 


THEORY 

It  is  often  difficult  to  determine  which  equation  will  yield  the 
best  results  when  solving  a particular  problem  using  an  energy  approach 
Therefore,  it  is  best  to  start  with  the  most  fundamental  equation  which 
encompasses  all  other  energy  approaches  to  solution.  Such  an  equation 
was  postulated  by  Sir  William  Rowan  Hamilton  and  was  called  by  him  the 
"Law  of  Varying  Action"  (Reference  9)*  It  can  be  mathematically  stated 
as , 

tl  ti 

6 f (T  + W)  dt  - il.  6q.  1 = 0 (1) 

where  T represents  the  total  kinetic  energy  of  the  system  and  W is  the 
total  work  of  the  forces  acting  on  or  within  the  system. 

1.  STRAIN-DISPLACEMENT  RELATIONSHIPS 

For  the  purpose  of  this  analysis,  the  plate  is  considered  to  be 
composed  of  an  isotropic  continuum  which  obeys  the  elastic  stress- 
strain  relationships.  In  addition,  the  assumptions  underlying  thin 
plate  theory  (plane  stress)  are  assumed  valid  for  the  vibrational 
analysis.  From  this  assumption  we  have, 
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which  reduces  the  stress-strain  relationships  to, 


o 
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<eX  + Vey> 


(1  - V2) 


(ey  + vex) 
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2(1  + V) 


(Yxy) 


Membrane  stresses  have  not  been  considered  for  this  analysis  (inexten- 
sional  plate  theory)  since  there  are  no  externally  applied  (thermal  or 
mechanical)  in-plane  loads.  From  the  assumption  of  small  deflections 
during  vibration,  the  corresponding  strain-displacement  equations 
reduce  to  a set  of  linear  equations  of  the  form, 
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3u 

9x 
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9v 

9y 
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Y » — + — 
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(3) 


With  the  two  assumptions,  plane  stress  and  small  deflections,  the 
in-plane  displacements,  u,  v,  can  be  expressed  as  functions  of  the 
out-of-plane  deflection,  w. 
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U = -Z 

3x 
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3w 

v = -z 

3y 


Substituting  Equations  4 into  Equations  3 and  combining  with 
Equations  2 yields. 
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(1 
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3y2 


xy 
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- (-  ^ ) 
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(5) 


2.  WORK 

For  the  conservative,  stationary  problem  of  the  free  vibration 
where  there  are  no  external  forces,  the  work  of  the  internal  stresses 
may  be  expressed  as  the  volume  integral  of  the  strain  energy  density 
function, 

W = - — / II  (a  e + a e + x y ) dx  dy  dz  (6) 
2 J JJ  * x yy  xy  'xy  7 


3 . KINETIC  ENERGY 

The  contribution  of  rotary  inertia  to  the  total  kinetic  energy 
has  been  found  to  be  insignificant  for  small  vibrations  of  a thin 
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plate.  Therefore,  the  kinetic  energy  can  be  expressed  as, 
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pQ  dx  dy  dz 


(7) 


As  a result  of  p,  w,  w satisfying  natural  continuity  conditions 
within  the  region  occupied  by  the  plate  material, 

3T  fi  rrr  • fi 

6q  = / / / pwdw  dxdydz 

1 't„  JJJ  ''o 


(8) 


Taking  the  variation  of  Equation  7 and  combining  with  Equation  8, 
as  in  Hamilton's  Law, 


L] " dt  ■ 5,1  L„  ~f!0  [Iff™ dx  dy  dz  ]■ 


dt  - 


Iff- 


(9) 


pw6w  dx  dy  dz 


Integrating  the  first  term  of  Equation  9 by  parts  with  respect  to 
time  and  carrying  out  the  necessary  algebraic  operations  yields. 
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pw6w  dx  dy  dz  • dt 


(10) 


4.  HAMILTON'S  ENERGY  EQUATION 

Substituting  Equations  6 and  10  into  Equation  1,  Hamilton's 
Law,  results  in  the  time  integral  of  the  principle  of  virtual  work, 
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ev  + o 6e  + t 6y  ) dxdydz 
x y y xy  ' xy  7 


dt  = 0 


The  resulting  equation  within  the  time  integral  in  terms  of  the  dis- 
placements of  the  system  under  consideration  is  well  known. 
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Assuming  simple  harmonic  motion  of  the  vibrating  plate  (free  vibra- 
tion), the  displacement,  w,  can  be  expressed  as, 


w(x,y,t)  = w(x,y)  sin  ut 


Substituting  this  expression  into  Equation  12  and  collecting  terms 


yields, 
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Recognizing  that  the  integral 


sin  utdt 


cannot  vanish  for  any  time,  tj_,  larger  than  tQ  , implies  that  the 
following  condition  must  exist. 
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dxdy  = 0 (15) 


Equation  15,  as  derived  is  applicable  to  only  simple  harmonic  motion 
of  the  system. 

5.  PLATE  COORDINATE  CONFIGURATION 

To  aid  in  satisfying  the  necessary  boundary  conditions,  the  plate 
has  been  modeled  as  two  discrete  elements,  one  of  which  embodies  the 
free  surface  of  the  crack.  The  coordinate  system  used  for  the 
analysis  of  the  combined  plate  elements  is  shown  in  Figure  1. 
Dimensional  characteristics  of  the  elements  served  to  model  the  root 
crack  configuration  for  a series  of  crack  depths  ( % chord) . The 
nondiraensionalization  performed  consisted  of  the  following: 
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nl  = n2  = y/a  ?i  = xl/bl  = x2^b2 
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AR-i  = a/bi 


Combined  Plate  AR 


ARo  = a/b. 


AR^  • AR2 


AR^  + AR2 


6.  DISPLACEMENT  FUNCTIONS 

Based  on  the  concepts  as  defined  by  Ritz  (Reference  2)  and  knowledge 
gained  from  References  4 through  8,  a set  of  admissable  functions  in  the 
form  of  a simple  truncated  power  series  served  as  the  deflection  function 
for  this  analysis.  The  functions  used  for  solving  the  crack  plate 
problem  are  expressed  in  the  nondimens ional  form. 


N M i i 

w1  UlS  n)  = r\)  Z l AtJ  n K1 

j=0i=0 


(16) 


N M 


irJ 


W (C2>  h ) = ^2’  n ^ 1 Bij  n ^2 

J=0i=0 


where  the  term  g(£,  n)  forces  satisfaction  of  the  prescribed 
geometrical  boundary  conditions.  The  most  general  form  of  this  func- 
tion for  the  coordinate  system,  Figure  1,  chosen  is.  Fig  1 


8<C.n>  - ?Nl  - O%01(1  - n)02 


(17) 
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The  specific  values  of  the  $s  and  0s  in  Equation  17  for  the  two 
element  solution  are. 

Element  1:  0^  » 2 <J>1  = 4>2  — = 0 

Element  2:  ■ $2  “ ®i  = = ® 


(18) 


For  the  line  of  commonality  between  the  two  elements  the  following 
continuity  conditions  are  required. 

Displacement:  (1,  q)  = W2  (0,  q) 


Chordwise  Slope:  (1,  n)  * W ^ (0,  n)  (19) 

x x 

Spanwise  Slope:  Wi  (1,  n)  = W (0,  n) 

V 2y 

The  author's  experience  (Reference  8)  in  performing  these  types  of 
calculations  for  the  case  of  beams  with  discontinuities  has  shown  that 
satisfying  only  the  slope  and  displacement  conditions  along  the  con- 
necting boundary  has  provided  excellent  convergence  on  the  higher  order 
derivatives  corresponding  to  the  moments,  shears,  and  forces. 

Satisfying  the  three  continuity  requirements,  Equations  19,  and  sub- 
stituting the  appropriate  "g"  functions.  Equation  18,  into  Equations 
16  yields. 


W 


1 


N M 

<{1’  n)  - E E n1  * h 
,1=0i=0 


J 
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N M 


w~  ($„,  n)  = e z a,  . n1  + 2 + aE  u.Jjln1'  + 


N M 


i + 2. 


2 2 


ij 


J=0i=0 


ij 


j=0i=0 


(20) 


N M 

S E B .yd 

j=2i=0  J 


Replacing  the  deflection  functions  in  Equations  18  with  Equations  20 
and  performing  the  necessary  operation  with  the  operator,  6,  where 


6W1  • E E 

JUok=0 


NMk  + ? N M , 

2 = S 5\4+aE  ^nk  + 2? 

•A  -0k=0  A=0k=0 


N M 
E Z 

■S  =2k=0 


2 \a  + 1 E 


results  in  the  set  of  equations  of  the  form. 


N M 

E E A, 
j=0i=0 


K 


+ B,,  (K2,,„,  - »2M2.jkil)J  ■ 0 


. . V IY_L  . . , - - A 1VLL  . , . * J T D \ a 

ij  L ijhj  lJkX  ij  ljkl 


(21) 


for  k = 

0, 

1. 

2 . 

. . M 

Jt  - 

0, 

1. 

2 . 

. . N 

3 = 

2, 

3 

. . N 

N M 


lK2TlJki  " * VK3W  - ‘Hja’  ] 


= 0 


(22) 


for  k = 0,  1,  2 
T = 2,  3 
^=2,3 


. M 
. N 
. N 


The  recursion  formulas  for  the  matrix  elements  are  presented  in  the 
Appendix.  Equations  21  and  22  can  be  expressed  in  matrix  form 


as , 


[Kl] 

[K2] 

- X2 

(MU 

(M3 

(1 

[K2]T 

(K3) 

(M2]  T 

m 

W 

= 0 (23) 


IK 


IM 


The  IK  and  IM  are  square  symmetric  matrices  of  order  M • (2N  - 2). 
For  a nontrivial  solution  to  this  set  of  homogeneous  equations,  the 
determinant  of  the  coefficient  matrix  must  vanish. 


[IK]  - X2  [1M]|  - 0 


(24) 
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SECTION  III 


L 


RESULTS 

Before  attempting  the  solution  to  the  nonuniform  cantilever  plate 
problem  using  the  method  of  Ritz,  a study  of  the  rate  and  type  of 
convergence  possible  was  performed.  The  uniform  cantilever  plate 
vibration  problem  served  to  demonstrate  the  convergence  possible 
using  this  approach.  A great  deal  of  data,  both  experimental  and 
analytical,  was  available  for  comparison.  Table  1 summarizes  the 
eigenvalue  comparisons  for  the  various  solution  schemes. 

A convergence  study  was  conducted  for  the  uncracked  configuration 
on  the  effects  of  varying  the  aspect  ratios  of  the  two  separate  elements 
while  maintaining  the  same  combined  aspect  ratio  of  2.  Table  2 
summarizes  the  results  of  this  study  for  four  different  element 
combinations.  As  depicted,  little  variation  was  noted  for  the  cases 
considered.  In  addition.  Table  3 shows  the  results  of  a study  con- 
ducted on  the  rate  of  eigenvalue  convergence  as  a function  of  the 
number  of  terms  used  in  the  truncated  power  series.  A study  of  this 
type  was  not  performed  on  the  effect  of  the  number  of  elements  in 
the  finite  element  solution  due  to  cost  restraints.  The  average  cost 
of  extracting  approximately  84  eigenvalues  and  their  corresponding 
eigenvectors  using  the  direct  approach  was  $15.00.  However,  the  cost 
of  extracting  the  first  10  eigenvalues  and  mode  shapes  using  the  finite 
element  program,  NASTRAN,  was  $150.00.  A study  was  performed  and  docu- 
mented (Reference  10)  which  showed  a sufficient  lower  bound  convergence 
for  a rectangular  cantilever  plate  with  an  aspect  ratio  of  2 using  the 
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same  number  and  type  of  finite  elements.  In  all  cases,  good  agreement 
was  obtained  for  the  uniform  cantilever  plate  solution  schemes,  Ritz, 


and  finite  element  when  compared  to  experimental  and  other  analytical 
solutions . 

With  the  good  correlation  between  the  experimental,  direct,  and 
finite  element  solutions  for  the  uniform  plate,  the  problem  of  the 
mixed  boundary  condition  plate  response  was  addressed.  The  crack 
was  introduced  into  the  analytical  solutions  as  a free  boundary 
with  no  surface  interactions.  For  the  experimental  test,  a narrow  cut 
was  made  along  the  root  of  the  plate  attachment.  The  plate  was  not 
removed  from  its  support  throughout  the  entire  testing  sequence.  An 
acoustical  siren  served  as  the  excitation  source.  Resonant  conditions 
were  defined  by  two  separate  methods. 

One  approach  employed  the  use  of  a laser  interferometry  technique 
to  visually  detect  the  resonant  mode  through  the  laser  light  inter- 
ference patterns.  The  second  method  involved  the  use  of  an  oscilloscope 
and  a piezo-electric  accelerometer.  Resonant  conditions  were  determined 
by  plotting  the  input  forcing  function  (sinusoidal)  against  the  output 
signal  as  relayed  by  the  accelerometer.  The  plate  used  for  obtaining 
the  experimental  data  consisted  of  an  aluminum  flat  plate  1/8"  x 3"  x 6" 
rigidly  constrained  at  one  edge  by  two  solid  steel  blocks.  These 
blocks  were  attached  to  an  air-damped  (floating)  table.  Poisson's 

3 

ratio  (v)  of  0.3,  mass  density  of  0.100  lb/in  , and  elastic  modules  (E) 
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7 2 

of  1.0  x 10  lbf/in  were  assumed  for  the  analytical  portion  of  the 
study. 

Figures  2 through  9 portray  the  eigenvalue  results  for  the  first 
eight  out-of-plane  vibration  modes.  The  data  depicts  the  response  as 
a function  of  root  crack  length  for  the  three  solution  schemes.  All 
eigenvalues  have  been  normalized  to  the  uniform  (uncracked)  configuration 
corresponding  to  each  solution.  In  all  cases,  the  experimental  frequen- 
cies fall  below  the  two  analytical  solutions.  Reasonably  good  agreement 
can  be  seen  in  most  instances  between  the  finite  element  and  energy 
solutions.  In  addition,  Figures  10  through  16  show  the  variation  in  mode 
shapes  as  a function  of  crack  length.  The  experimental  mode  shapes  con- 
sist of  the  holograms  taken  from  the  laser  setup  previously  mentioned. 

The  finite  element  mode  shape  plots  are  contour  surface  diagrams,  whereas 
the  energy  solution  generated  plots  represent  positive  and  negative 
out-of-plane  displacements. 
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SECTION  IV 


CONCLUSIONS 

The  results  of  this  work  have  demonstrated  that  a direct  and 
economical  solution  to  the  nonuniform  plate  vibration  can  be  obtained. 
Excellent  agreement  was  shown  for  the  uniform  plate  vibration  solution 
when  compared  to  those  obtained  from  literature.  Comparison  of  the 
results  for  the  nonuniform  vibration  problem  has  indicated  some 
differences  between  the  theoretical  and  experimental  solutions. 

Similar  results  were  obtained  for  both  the  finite  element  and  energy 
solutions.  Figure  17  portrays  the  relative  comparison  between  the 
uniform  and  nonuniform  deviations  from  the  experimental  solution. 

As  shown,  a marked  increase  in  the  percent  difference  occurs  for  the 
cracked  plate  solution.  This  difference  may  possibly  be  explained  by 
either  a breakdown  in  the  assumptions  associated  with  thin  plate 
vibration  theory  or  the  approach  used  in  obtaining  the  experimental 
results.  Both  results  indicate  the  same  type  of  frequency  decay  with 
increasing  crack  length,  but  differ  slightly  in  the  absolute  value  of 
this  decay.  Additional  experimental  testing  and  analysis  are  necessary 
to  resolve  the  problem  of  the  noted  difference. 

This  paper  has  presented  an  essential  step  towards  the  under- 
standing of  the  effect  of  cracks  on  quasi-blade  response.  It  is 
apparent  from  the  results  that  further  work  in  the  area  of  nonuniform 
blade  response  is  warranted  before  consideration  is  given  to  extending 
the  damage  tolerance  criteria  to  gas  turbine  engine  airfoils. 


( 
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ROOT  CRACK  LENGTH  ( % CHORD  ) 


FIGURE  10.  ROOT  CRACK  MODE  SHAPE:  MODE  1 


FIGURE  13.  ROOT  CRACK  MODE  SHAPE:  MODE  4 


FICURE  15.  ROOT  CRACK  MODE  SHAPES:  MODE  6 


ROOT  CRACK  LENGTH  ( % CHORD 


FIGURE  16.  ROOT  CR.<\CK  MODE  SHAPES:  MODE 
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FIGURE  17,  EIGENVALUE  RATIO  DEVIATIONS  FROM  EXPERIMENTAL  SOLUTION 
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I12ijkA 


(i  + 2) (i  + 1) (A) 

(i  + k + 3) 
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ijkA 


(j)(i  + 2)  (k  + 2) 
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I14iJkJt 
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I15ijkA 
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(i  + k + 3)(j  + A - 1) 


I17ijkA 
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(i  + k + 3)(j  + A - 1) 
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APPENDIX  (CONTINUED) 
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ijkl 
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(i  + k + 5)(j  + X + 1) 

1 

(i  + k + 5) 
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j-i 

(i  + k + 5)  (3) 

1 

(i  + k + 5)  ( A + 1) 
j 

(i  + k + 5) (1  + 2) 

1 
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APPENDIX  (CONTINUED) 
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